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We propose a feasible route to engineer one and two dimensional time reversal invariant (TRI) 
topological superconductors (SC) via proximity effects between nodeless s± wave iron-based SC and 
semiconductors with large Rashba spin-orbit interactions. At the boundary of a TRI topological 
SC, there emerges a Kramers pair of Majorana edge (bound) states. For a Josephson 7r-junction we 
predict a Majorana quartet that is protected by mirror symmetry and leads to a mirror fractional 
Josephson effect. We analyze the evolution of Majorana pair in Zeeman fields, as the SC undergoes a 
symmetry class change as well as topological phase transitions, providing an experimental signature 
in tunneling spectroscopy. We briefly discuss the realization of this mechanism in candidate materials 
and the possibility of using s and d wave SC and weak topological insulators. 
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Introduction. — Broken symmetry and topological or- 
der are two fundamental themes of condensed matter 
physics. The search for topological superconductors (SC) 
[1-3] is fascinating, as gauge symmetries are sponta- 
neously broken in the bulk and gapless Andreev bound 
states (ABS) can be topologically protected at order pa- 
rameter defects, hosting Majorana fermions. Majoranas 
are immune to local noise by virtue of their nonlocal topo- 
logical nature and thus give hope for fault-tolerant quan- 
tum computing [4]. The rise of topological supercon- 
ductivity (SC) has been expedited by recent proposals 
[5-9] that hybridize ordinary SC's with helical materials, 
with the help of magnetic perturbations. Using proxim- 
ity effects, electrons in a single helical band at the Fermi 
energy form conventional Cooper pairs, whose conden- 
sation realizes a spinless chiral p wave SC in its weak 
pairing regime, i.e., a topological SC with broken time 
reversal symmetry (class D). Unique signatures, includ- 
ing zero bias conductance peaks, anomalous Fraunhofer 
patterns, and fractional Josephson effects, are starting to 
be observed in these systems [10-14]. A completely dis- 
tinct family (class Dili) of time reversal invariant (TRI) 
topological SC's was proposed based on a mathematical 
classification of Bogoliubov-de Gennes (BdG) Hamiltoni- 
ans [15-20]. Cu x Bi2Sc3 [20] and Rashba bilayers [21] are 
possible candidates, if a fully gapped inter-orbital spin- 
triplet pairing dominates. This seems challenging since 
exotic interactions are required. Experimentally [22-27], 
the nature of the pairing remains very controversial. 

A more ambitious goal is to realize TRI topological SC 
without exotic electron-electron interactions and without 
magnetic fields. Here we propose a feasible route to uti- 
lize proximity effect devices which combine Rashba semi- 
conductors (RS) and nodeless iron-based SC's. Below its 
transition temperature the SC provides the RS a s± wave 
spin-singlet pairing potential that switches sign between 
the T and M points [28-30]. TRI topological SC is re- 
alized when the chemical potential is adjusted to make 
the inner and outer Fermi surfaces feel pairing poten- 
tials with opposite signs. At a boundary of the 2D (ID) 




FIG. 1. Left panel: sketch of the proximity devices pro- 
posed in the text. A Majorana Kramers pair emerges at 
the boundary when the TRI SC becomes topological. Right 
panel: sketch of the Josephson junction described in the text. 

TRI topological SC, a Kramers pair of Majorana edge 
(bound) states emerge as localized midgap states. For a 
Josephson 7r-junction a Majorana quartet is protected by 
mirror symmetry, leading to a mirror fractional Joseph- 
son effect. We analyze how the Majorana pair evolves 
in Zeeman fields, as the SC undergoes a symmetry class 
change and topological phase transitions, providing an 
experimental signature in tunneling spectroscopy. 

2D TRI topological SC. — We first introduce a minimal 
model on a square lattice to characterize 2D TRI SC: 

H = -t c t* c jc ~ iX n J2 C L( CT a/3 x dij) z C j0 

<ij>,a <ij> 

+ Ao £(i c l + h-c) + A x £ ( C y u + h.c.) .(1) 

i <ij> 

Here t is the nearest neighbor hopping and er are the 
Pauli matrices of electron spin. The second term arises 
from the Rashba spin-orbit interactions. Note that dij is 
a unit vector pointing from site j to site i and we assume 
Ar, > 0. Ao and Ai, induced by the proximity effect, 
lead to a combined s± wave pairing potential. It is more 
convenient to write the BdG Hamiltonian: 

^fc d ° = [-2t(cos k x + cos k y ) + h% -fj]r z + A k t x 
ft.}? = 2Ar,(sui k x a~y — sin k y o~ x ) 
Afe = A + 2Ai (cos k x + cos k y ) , (2) 
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where \i is the chemical potential and r are the Pauli ma- 
trices in Nambu particle-hole notation. is a s± wave 
singlet pairing potential that switches signs between the 
zone center T (0, 0) and the zone corner M (tt, it) when 
the order parameters satisfy |A | < 4|Ai|. As we show 
in Fig. 1(a) and note below, A& could be provided by a 
nodeless iron-based SC on which the Rashba layer is de- 
posited. For convenience we have chosen a gauge to pin 
the overall phase of s± wave SC at zero and assumed 
< |Ao| < 4Ai hereafter. "H^ dG has time reversal 
(O = —iayK) and intrinsic particle- hole (S = a y T y K) 
symmetries. We obtain the energy dispersion 

E% dG = ± yj [2i(cos k x + cos k y ) + ii ± e\\ 2 + A 2 , (3) 

where = iX^^Jsm 2 k x + sin 2 k y is the Rashba en- 
ergy. Afc has a closed nodal line, i.e., cosk x + coak y = 
—A /(2A 1 ), in the first Brillouin zone. At the nodal line, 
£BdG = ± ( M _ £o ± e R) with eo = tAp/Ai, and e£ has 

the maxima e^ ox = 2Ar-\/2 — A 2 ,/ (8A 2 ) and the minima 
4 <n = 2ArV|A /A 1 |-A 2 1 /(4A 2 ). 

In both 2D and ID, the Z2 topological invariant [15- 
18] of a TRI SC is determined by whether the pairing 
potential has a negative sign on odd number of Fermi 
surfaces each of which encloses a TRI momentum [17]. 
As shown in Fig. 2 and summarized in Table I [31], the 
phase of the hybrid SC depends on the chemical potential 
M- For the case of e* in < |/x-e | < e* QX , H^ dG describes 
a nodal SC. When \p — e | > e„ al , the SC is fully gapped 
but in the trivial (y — 0) phase since A& has the same 
sign on both Fermi circles. When — e 1 < eJJ in is satis- 
fied, the pairing potential switches sign between the two 
Fermi circles [32] , and consequently the hybrid system re- 
alizes a TRI topological SC (y — 1). The energy window 
for tuning the system into the v = 1 state has the size 
of 2eJ^ in with an optimized value 4Ar, at A = ±2Ai. 
For the v = 1 state helical Majorana edge states emerge 
at the boundary, as shown in Fig. 2, the spectrum of a 
infinite ribbon described by Eq. (1). The Majoranas at 
k = 7T (0) for sgn(A /Ai) = 1 (—1) are protected by 
time reversal and particle-hole symmetries. 

ID TRI topological SC. — By turning off all the k y 
terms Eq. (2) models a ID Rashba nanowire deposited 
on a nodeless s± wave SC. When the two s± wave or- 
der parameters satisfy |Aq| < 2Ai, the pairing potential 



TABLE I. Summary of the Z2 classification of the hybrid TRI 
SC in class DHL eo, eJJ,^, and e^„ are defined in the text. 



Phase 


Two Dimension 


One Dimension 


v = 1 




[/Lt - e < fm 


Nodal 


^min — IM ^0 1 — ^-max 


j/i — eo = ef n 


v = 


\^ ^0 1 ^ ^max 


\fJ, — e | > ef n 




FIG. 2. Upper panels: the two Fermi surfaces (blue and 
green) of the single-particle bands for the v — 1 state; the 
closed nodal line (red) of Afe, separating two regions in which 
Afe has opposite signs. Lower panel: BdG spectrum of a 2D 
ribbon as a function of k for the v = 1 state with p, = eo. 
The red and green lines indicate the helical Majorana edge 
states. We choose parameter values: t — 10, Ar = 5, and 
|A | = Ai = 2. (a) A > and (b) A < 0. 



switches sign between the two TRI momenta and ir. In 
ID, the closed nodal line of A;, is shrunk to two nodes at 
k = ± arccos(— Ao/2Ai). At the nodes, the Rashba en- 
ergy is e^; = 2Ar\/i — Aq/(4A 2 ), and thus a proximity 
induced ID TRI nodal SC is identified for \i = e ± e^- 
When \fi — e | < ej^, a positive pairing is induced for the 
inner pair of Fermi points while a negative pairing for 
the outer pair, realizing a ID TRI topological SC. In the 
case of — eo I > e^, the hybrid system becomes a trivial 
SC that is adiabatically connected to the vacuum state. 

At each end of a ID TRI topological SC, there emerges 
a Kramers pair of Majorana bound states (MBS). With- 
out loss of generality, in the rest of this paper we will 
set Ao = for the ID case and thus < 2Ar is the 
criterion for the v = 1 state, as shown in Fig. 3(a). The 
cyan line denotes four degenerate MBS's independent of 
Aj. Further investigation of their wavefunctions shows 
that these four MBS's form two Kramers pairs localized 
at the opposite ends of nanowire. This verifies our ana- 
lytical results summarized in Table I. 

Mirror Fractional Josephson effect. — Consider the lin- 
ear Josephson junction in Fig. 1(b), in which a Rashba 
nanowire is deposited on a larger s± wave SC ring, and 
the phase difference (j> — (2e/h)$ across the junction 
is controlled by the magnetic flux $ through the ring. 
Fig. 3(b) shows the spectrum of ABS's as a function of </> 
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FIG. 3. (a) BdG spectrum of the ID TRI SC as a function 
of /i. (b) Spectrum of ABS's in the junction as a function 
of <j>. The cyan lines are four-fold degenerate, denoting two 
Kramers pairs of MBS's at opposite ends. The red and green 
lines are doubly degenerate, denoting two pairs of ABS's in 
the junction. We choose parameter values: t = 10, Ar = 5, 
Ai = 2, and Ao = 0. fj, = is used in (b). 

when the physical separation between the ends of SC ring 
is small. The four-fold cyan line shows the appearance 
of a pair of MBS's at each end. The red and green lines, 
both doubly degenerate, represent two pairs of ABS's in 
the junction. When cf> ^ tt, a finite \x lifts the degener- 
acy. However, protected by time reversal, particle-hole 
and mirror (M y = —io~ y ) symmetries, their crossing at 
E = when <f> = tt is a four-fold degeneracy of special 
significance. This 7r-j unction is in sharp contrast with 
its Z2 counterpart in the class without mirror symme- 
try [8, 18, 34]. Splitting the Majorana quartet at <f> = tt 
would break mirror and/or time reversal symmetry. To 
understand this topological twist in Fig. 3(b), we lin- 
earize Eq. (2) near the Fermi energy: 

Heff = (vk x a y - c) st z ~ fir z + Use , (4) 

Use = A s [cos - t x + sgn(x) sin ^ r y ] , (5) 

where s — ± denote the inner and outer bands with 
opposite spin helicities and c (0 < c < /1) lifts their 
degeneracy [33]. The mirror symmetry allows one to 
label the bands with the eigenvalues of M y , a = ±i. 
We find the ABS dispersions e a (</>) = Acos(</>/2). Note 
that the perfect normal state transmission and the in- 
dependence on fi, c and a are artifacts of the simpli- 
fied model. We can define Bogoliubov operators T a ± 
that satisfy T a = T a+ = r^_ because of the particle- 
hole symmetry. The low-energy Hamiltonian is thus 

H = E a £ a(^)( r l r a - I) = 2i£a e "(^)7a??a where 
7« = (It + r Q )/2 and ^ = i(r], - T a )/2 are the Ma- 
jorana operators. For each band, N° = rj,r a = 0, 1 
distinguishes two states and coupling them requires an 
process that changes N®. Due to the Cooper pairing, 
the total charge is not conserved. However, the fermion 
parity of each band N a mod 2 is conserved, as the mir- 
ror symmetry does not allow scattering between the two 



bands. This mirror fermion parity conservation forbids 
the mixing among the four ABS's in the junction and 
therefore protects their crossing at zero energy. cf> acts 
like a defect and parameterizes the mirror fermion parity 
pump. Although Eq. (4) is invariant under <5$ = h/2e, 
the global Hamiltonian is physically distinct. When a 
flux h/2e is threaded through the SC ring, <fi is advanced 
by 27r, 7 Q — > 7 Q while r\ a — > —r] a , and a unit of fermion 
parity is transferred between the two bands resolving the 
fermion parity anomaly for an individual band. 

In response to the change of phase, the populated 
ABS's carry supercurrents I a ± = ±I a through the de- 
vice, whereas the states in the continuum has negli- 
gible contributions. For the v = 1 case, we obtain 
I a = (2e/h)de a /d(f> ~ (— eA/h) sm(cf>/2), which is maxi- 
mized at <j) = 7r in sharp contrast to the v — (conven- 
tional) case. In the absence of mirror symmetry break- 
ing, there is no transition among I a ±, signaling a mirror 
fractional Josephson effect with 47r periodicity. 

Evolution of Majorana pair in Zeeman field. — When 
one helical band is removed from the Fermi energy, a 
Rashba nanowire proximity coupled to a s wave SC is a 
topological SC with broken time reversal symmetry, sup- 
porting a single MBS at each end. Realizing such a topo- 
logical phase requires the Zeeman field V z o~ z (or V z a x ) 
and the chemical potential /x to satisfy: 4A^ + (|/i|— 2t) 2 < 
V z < 4Af + (\fi\ +2t) 2 . Since only one helical band is 
present at the Fermi energy this also occurs even if the SC 
is s± wave as long as the hybrid SC remains fully gapped 
(/i 2 7^ 4A^ + V z ). The latter condition is guaranteed for 
the TRI v = 1 state as it satisfies < 2Ar. It is thus 
intriguing to investigate how a MBS pair evolves in the 
Zeeman field, as the bulk SC undergoes a symmetry class 
change and topological phase transitions. Fig. 4(a) shows 
the evolution in the case of \i ^ 0. When V z is turned on, 
without gap closing, the topological SC in the TRI class 
becomes a trivial SC in the class without time reversal 
symmetry. Two topological phase transitions occur at 
V 2 = 4A 2 + (|/x| ± 2i) 2 where the gap closes. The v = 1 
state in the new class is realized between the transitions. 
At one end, as V z is tuned up, one MBS disappears at the 
first transition while the other persists in the v = 1 state 
and enters the bulk continuum at the second transition. 
For the special case [i = 0, the two transitions merge into 
one and the new v = 1 state does not appear. 

As implied by the zero energy black lines in Fig. 4(a), 
it seems baffling that the Majorana pair is robust against 
the Zeeman field before the first transition occurs. We 
emphasize that these Majoranas are not topologically 
protected because the bulk SC is trivial in the new class. 
However, their robustness can be understood using the 
effective model described by Eq. (4) with 

Use = A s Tx 6{a - \x\) + At x 6(\x\ - a) , (6) 

where x = ±a are the locations of the two ends of SC 
and A — > 00 reflects the infinity mass of vacuum. Eq. (6) 
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FIG. 4. Evolution of a MBS pair in Zeeman fields, (a) 
In a a z (cr x ) field; (b) in a a y field. The cyan dots (lines) 
have two-fold (no) degeneracy, indicating the appearance of a 
MBS pair (single MBS) at the end of SC. The black lines at 
zero energy are doubly degenerate. The red and green lines 
indicate the Zeeman splitting of two zero energy states with 
opposite fcrmion parity. We choose parameter values: y, = — 5 
in (a), fi — in (b), and others are the same as in Fig. 3. 



incorporates the correct boundary condition [35] of TRI 
topological SC. Solving the boundary problem, we find 
that a pair of MBS's at one end forms two states with 
a y = T y = ±1 while the other pair at the opposite end 
forms two states with a y = —r y = ±1. Clearly, a o~ x or 
o~ z Zeeman field only couples states with opposite o~ y fla- 
vors at opposite ends, whose wavefunction overlap decays 
exponentially on the SC length. Consequently, none of 
the Majoranas can be passivated away from zero energy 
until the first topological phase transition. However, as 
shown in Fig. 4(b), a small a y field can Zeeman split the 
two zero energy states with opposite fermion parity at 
each end, revealing the topological triviality of bulk SC. 
The a y field also closes the SC gap. 

The evolution of a Majorana pair in a Zeeman field 
provides a smoking gun for the identification of a TRI 
topological SC in tunneling spectroscopy. The Majorana 
pair induces resonant Andreev reflection producing a zero 
bias conductance peak [36-38] of 4e 2 jh. This peak per- 
sists for small a x and a z fields, reduces to 2e 2 /h after 
the first topological quantum phase transition, and dis- 
appears at the second one. In contrast, a small o~ y field 
not only Zeeman splits but also reduces the peak. 

Discussions. — Unlike a d xy wave SC, where it is im- 
possible to induce a full gap at the Fermi surface centered 



at a TRI momentum, a nodeless s± or d x 2_ y 2 wave SC 
allows the pairing potentials on Fermi surfaces centered 
at (0, 0) and (tt, it) or at (0, n) and (n, 0) to have opposite 
signs. Therefore, it would be interesting to consider 2D 
TRI topological SC in a proximity device that hybridize 
a weak topological insulator with two surface states and 
a nodeless s± or d x 2_ y 2 wave SC [31]. Although the ro- 
tational symmetry of RS makes the hybrid 2D d x 2_ y 2 
wave SC nodal, it is still possible to use them to engineer 



a ID TRI topological SC. This has been recently sug- 
gested by Wong and Law [39], however, the d x 2_ y i wave 
heavy fermion SC they chose is nodal, which provides a 
decay channel for MBS's and destroys their topological 
robustness. One may also wonder whether it is possible 
to build a hybrid TRI topological SC using a s wave SC 
if the proximity effect is strongly orbital-dependent. Un- 
fortunately, this seems not plausible for the reason that 
we elaborate in the supplementary materials [31]. 

Our proposal is experimentally feasible, since the nec- 
essary ingredients and the required technologies are all 
well established. It has been widely accepted that, at 
least for iron pnictides, there are a large family of node- 
less s± wave SC's [28-30], though their pairing mecha- 
nism is still under lively debate. An iron-pnictide with 
closer electron and hole Fermi surfaces is preferred. Ao,i 
can be adjusted by doping or changing materials. Iron- 
pnictides consist of square layers that are stacked in 
tetragonal structures and it is better to choose RS with 
cubic or tetragonal crystal structures, e.g., Au, Ag and 
Pb, in order to match their Brillouin zones (BZ). A small 
lattice incommensurability may blur and thus effectively 
broaden the pairing potential in the BZ of RS, which may 
even be helpful as long as the SC remains nodeless. 

Our proposal for a realization of TRI topological SC's 
in ID and 2D constitutes four critical advances: (i) there 
is no need for magnetic perturbations or exotic interac- 
tions, simplifying the experimental setup; (ii) the SC gap 
can reach more than 15 meV, raising the critical temper- 
ature of topological SC; (hi) it is likely that the presence 
of time reversal symmetry mitigate the role of disorder 
via Anderson's Theorem, making the topological SC even 
more robust; and (iv) irrelevant bands [40] are completely 
absent, allowing large tunability in feasible materials. 
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I. 2D time reversal invariant topological Superconductivity 

We further illustrate the 2D Z2 classification described in the main text and summarized in Table I. The energy 
dispersion at the nodal line of pairing potential reads 

E™ G = ± („ - e ± e*) , (1) 

where is the energy associated with the Rashba spin-orbit couplings. At the nodal line, has the maxima 
e* ax = 2A RV /2 - A2/(8A2) and the minima ef nm = 2A RV /|A /Ai| - Ag/(4A?). The lower panel of Fig. 1 depicts 
the energy dispersions of the two single-particle bands at the nodal line of pairing potential, (i) For the case of 
e min — I M — e 1 < e^ ax , "Hf dG describes a nodal superconductor, as Eq. (1) can be zero for some momenta, (ii) When 
e.o\ > e^ aa ., the superconductor is fully gapped but in the trivial (y = 0) phase, because A& has the same sign at 
both Fermi circles, (iii) Only when \[i — e | < is satisfied, the pairing potential switches sign between the two 
Fermi circles, and consequently the hybrid system realizes the time reversal invariant topological superconductivity 
(i/ = l). 




k x ■ n k x 7i 



FIG. 1: Upper panels: the two Fermi surfaces (blue and green) of the single-particle bands for the v = 1 state; the closed 
nodal line (red) of A^, separating two regions in which A*, has opposite signs. Lower panels: the energy dispersions of the two 
single-particle bands at the nodal line of Afe in the region < k x > n/2 and k y > 0, where eo is defined as the energy origin; /j, 
lies in yellow, cyan and white regions correspond to v — 1, v = and nodal states in Table I. 
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II. Using a weak topological insulator 

A nodeless s± or d x 2_ y z wave superconductor allows the pairing potentials on Fermi surfaces centered at (0, 0) and 
(-7T, 7r) or at (0, 7r) and (tt, 0) to have opposite signs. Therefore, it is promising to build 2D time reversal invariant 
topological superconductivity in a proximity device that hybridize a weak topological insulator with two surface states 
and a nodeless s± or d x 2_ y 2 wave superconductor. This two possibilities are sketched in Fig. 2. 




-7T 71 -7T 7T 

FIG. 2: In each panel, the blue and green circles denote two Fermi surfaces of the surface state of a weak topological insulator, 
and the red line separates two regions in which A*, has opposite signs. 



III. Possibility of using a s wave superconductor 

To elaborate this point, we use Eq. (4) with c = and s replaced by s z . s z — ±1 denote the top/bottom surface 
state (or layer) in the case of strong topological insulators (ref. 35) or bilayer Rashba films (ref. 24). ms x r z can be 
understood as a small trivial gap at k = due to the tunneling between the (top-or-bottom) orbital-pseudospins. In 
the weak tunneling limit, the proximity induced pairing potential has the form of 

Tt x = [(4 + t%)/2 + (4 - t%)/2 s z + t T t B s x ] t x , (2) 

where tj- (tg) is the tunneling amplitude between the superconductor and the top (bottom) orbital. We have assumed 
tx and ts are real without loss of generality, and will use Ti tZtX to denote the three coefficients quadratic in t in the 
bracket. Note that only the s z t x pairing, breaking inversion symmetry, gives rise to the topological superconducting 
state in class DHL However, such a state will compete with the trivial superconducting state. Focusing on %sc, 
the superconducting gap is proportional to \Tj — \fT% + T 2 | which vanishes. Therefore the hybrid superconductor is 
nodal. Even if the inter-orbital pairing is not allowed, i.e., T x = 0, Tj > \T Z \ is generally valid and therefore pins the 
hybrid superconductor to the trivial state. When the s wave superconductor is integrated out in the standard second 
order perturbation theory, the renormalized pairing potential reads 

^ = VA 2 - u? ^ + VA 2 - ^ ^ 

with r^v = itNqVT. Here V and No are the tunneling volume and the density of states of the superconductor in the 
normal state at the Fermi energy. Interestingly and obviously, the conclusions are unchanged at the Fermi energy in 
this more rigorous approach. Similar analysis and conclusions can be generalized to the case with both orbital and 
spin flavors, in which the spin-triplet pairing o~ z s y T x leads to a topological state. 
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